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Abstract 


Introduction 


Mean velocity and turbulence measurements were 
conducted on the three-dimensional fountain flow- 
field generated by the impingement of two axisym- 
metric jets on a ground plane with application to 
vert i cal- take-off and landing (VTOL) aircraft. The 
basic instantaneous velocity data were obtained 
using a two-component laser Doppler velocimeter in 
a plane connecting the nozzle centerlines at dif- 
ferent heights above the ground emphasizing the jet 
impingement region and the fountain upwash region 
formed by the collision of the wall jets. The 
distributions of mean velocity components and 
turbulence quantities, including the turbulence 
intensity and the Reynolds shear stress, were 
derived from the basic velocity data. Detailed 
studies of the characterist i cs of the fountain 
revealed self-similarity in the mean velocity and 
turbulence profiles across the fountain. The 
spread and mean velocity decay character i sties of 
the fountain were established. Turbulence inten- 
sities of the order of 50^ were observed in the 
fountain . 


Nomenclature 

D Exit internal diameter of the nozzle or 

the exit jet diameter 

H Height of the nozzle exit above the ground 

plate 

Re Reynolds number based on the exit jet 

diameter and the exit jet velocity 
S Center distance between the jets 

U Mean velocity in the streamwise direction 

V Mean velocity in the cross-stream 

direct ion 

u Fluctuating component of the mean velocity 

in the streamwise direction 

v Fluctuating component of the mean velocity 

in the cross-stream direction 

Uj Jet centerline velocity at the exit 

U The local maximum streamwise velocity in 

m 3 x 

the fountain 

X, . The fountain half width where U = -Ir U 

v 1 /2 _ • 4 . ......... 2 . max 

X Distance parallel to the line connecting 

the nozzle centerlines 

Y Distance perpendicular to the line con- 
necting the nozzle centerlines and paral- 
lel to the ground plane 

Z Perpendicular distance from the ground 

plane 


The impinging lift jets of a vertical take-off 
and landing (VTOL) aircraft hovering in ground 
proximity produce complicated three-di;. ^nsional 
flowfields. These 3 - D flowfields shown schemati- 
cally in Fig. 1 for the case of two-jets, involve 
strong interactions between the impinging lift-jet 
streams, the airframe surface and the ground. An 
important feature of this flowfield is the fountain 
upwash flow generated by the colliding wall jets. 
The fountain is fan-shaped (Fig. 1), spreading 
radially in all directions with increasing width 
away from the ground. The impingement of the 
fountain on the aircraft increases lift, elevates 
skin temperatures and causes possible reingestion 
into the inlets. Detailed studies of the structure 
and development of the fountain upwash flow are 
thus essential for the successful design of effi- 
cient VTOL aircraft. Toward this objective, a 
unique set of velocity and turbulence data on an 
axisymmetric twin-jet fountain flow are presented. 


Background 

References 1-4 report some of the earlier 
measurements of fountain upwash flow properties. 
Reference 1 reports mean velocity measurements in 



1 . Lift jet flow 

2. Jet impingement region 

3. Wall jet flow 

4. Fountain formation region 

5. Fountain up-wash flow 

6. Wall jet interaction stagnation line 

7. Entrainment 

8. Ground plane 

9. Blocking surface gp 2 i-ii 2 *-j 


Fig. 1 Schematic illustration of twin-jet 
impingement flow. 
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the wall jet and the fountain for twin axisymmetric 
jet impingement using hot-wire techniques for S/D * 
12 and H/D = 5. The fountain formation and devel- 
opment characteristics were obtained for a range of 
the important parameters that influence the flow. 

The flow angularity in the fountain in the plane 
connecting the nozzle centerlines was reported to 
vary randomly. Reference 2 presents pitot pressure 
profiles in the fountain for twin-axisymmetric jet 
impingement with S/D ■ 6 and 2 < H/D £ 6. Strong 
fluctuations were reported in the upwash flow 
properties, and the magnitude of these fluctuations 
increased with jet height above the ground. Also, 
highly fluctuating flow angularity, with frequent 
complete velocity reversals, was reported in the 
regions outside of the main upwash fan and the 
shear layers of the incident jets, and sometimes 
even in the central upwash region. Therefore, it 
was recommended that pitot probes be used in the 
measurement of mean flow properties in the fountain 
upwash . 

Reference 3 reported measurements in a two- 
dimensional fountain upwash, indicating high tur- 
bulence levels and spreading rates in the fountain. 
The first set of careful measurements conducted in 
a two-dimensional fountain upwash using x— y i re hot 
film anemometer were presented by Gilbert, who 
reported detailed mean velocity and turbulence 
measurements, including decay and spread charac- 
teristics of the fountain. The observed levels of 
turbulent intensities in the fountain were similar 
to those in an ordinary two-dimensional free jet; 
however, high growth rates were observed in the 
fountain. The fountains in References 3 and 4 were 
generated by two isolated, two-dimensional opposing 
wall jets, thus eliminating the jet impingement 
region as part of the fountain formation process. 

Reference 5 presents mean velocity and tur- 
bulence measurements using hot-wire techniques in a 
fountain generated by the impingement of two-axi- 
symmetric jets, with no definite conclusions on the 
fountain turbulence structure. References 6 and 7 
describe fountain behavior for close nozzle spacing 
and also indicate the effect of the blockage by the 
probe support. Measurements reported in Refs. 8-11 
do not reveal any additional features of fountain 
flows . 

Although several investigations of the fountain 
flow have been carried out, interpretations of the 
measurements vary widely, primarily because of the 
difficulty in measuring highly unsteady flows using 
hot-film and pitot-probe techniques. Computational 
codes require a better definition of the fountain 
flow and its turbulence structure to make reliable 
predictions. Therefore , a need existed for a 
reliable, detailed mean-velocity and turbulence 
data base in a realistic 3~D fountain generated by 
the impinging jets using a suitable measurement 
techni que . 

Measurement Technique and the Working Medium 

The measurement technique for studying the 
fountain flow was required to meet the following 
specifications: a) be nonintrusive so that no flow 

disturbance is caused, especially in the case of 
closely spaced jets, b) be able to sense the direc- 
tion of flow in the recirculat ing regions, c) have 
a linear response, essential for accurate measure- 
ments in areas where the turbulence intensities are 


high such as the fountain. A laser Dopp^r^ velo- 
cimeter ( LDV ) meets these specifications * ; in 

addition, LDV responds to a specific velocity 
component in 3“D flows and measures velocity di- 
rectly without the need to correct for temperature 
effects . 

As a working medium, water offers specific 
advantages compared to air; the tracers suitable 
for flow visualization in water are more numerous 
with better light-reflecting characteristics, and 
aerodynamic phenomena can be observed at a rela- 
tively slow speed for the same Reynolds number and 
model scale because of the difference in kinematic 
viscosities of air and water. Moreover, for laser 
Doppler velocimetry, the seeding of the flow is 
normally unnecessary in water because the natural 
suspended particles act as light scatterers. If 
seeding is necessary, it can be accomplished easily 
in water than in air by adding neutral-density 
plastic particles of the proper size. 

Experimental Apparatus 

Figure 2 shows the jet-impingement facility 
used for conducting laser-Doppler-velocimeter (LDV) 
measurements with water as the working medium. The 



1 Main plexiglass tank 

2 Header tank 

3 Nozzle units 

4 Ground plate 

5 Turbine flow meter 

6 Motor-operated 
flow- control valves 

7 Shut-off valves 

8 Settling chamber of 3 

9 Control pand 

10 Traversing mechanism 
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Fig. 2 Jet Impingement facility. 


ORIGINAL PACE »s 
OF POOR QUALITY 


OF POOR QUALITY 


primary components are (1) the main Plexiglas tank, 
(2) the header tank, (3) nozzle units, and (4) the 
ground plate. The large transparent main tank 
(1 .88 x 1 .52 x 1.37 m) has an approximate capacity 
of 3800 liters. The pressurized stainless steel 
cylindrical header tank, 91 . 4-cm diam., supplies 
water to the nozzle units and acts as a plenum to 
damp fluctuations produced by the pump. 

The function of the Plexiglas nozzle units 
(Fig. 3) is to reduce turbulence and mean-velocity 
nonuniformit ies in the existing jet-flow to accept- 
able levels. Each nozzle unit consists of a flow 
distributor (diffuser) for decelerating the flow, a 
honeycomb and screens to establish a uniform flow 
with low turbulence, and a nozzle (2.54-cm exit 
diam., 16:1 contraction ratio) for accelerating the 
flow. The ground plate is held at a fixed height 
above the bottom of the main tank and acts as an 
impingement surface for the jets with sufficient 
edge clearance for passage of the flow. The posi- 
tion of the nozzle units is adjusted by a travers- 
ing unit that can accommodate as many as four 
nozzle units and mounted on the main tank. 

The circulation system draws off water from 
under the ground plate and pumps it into the header 
tank, thus supplying water to the jets through a 
series of flow-control devices. With the available 
pumping capacity, jet Reynolds numbers up to 
200 000 can be obtained in a typical twin-jet 
(2.54-cm-exit diam.) impingement configuration. 

The three-dimensional fountain flow is visualized 
through a series of two-dimensional images. Fluor- 
escein-sodium, a fluorescent dye, is injected into 
the jet flow, which fluoresces bright yellowish- 
green when the desired cross-section is illuminated 
with a thin (1 mm) sheet of 488 nm wavelength light 
from an Argon-ion laser. 

A two color (two-component) TSI laser Doppler 
velocimeter system (Fig. 4) is used in the dual- 
beam off-axis backward scattering mode with Bragg- 
cell frequency shifting. The probe volume is 
positioned at the required location using a remote- 
ly driven X-Y-Z traversing unit. The TSI counter- 



type signal processors were used to convert the 
Doppler signal into a form suitable for recording 
on a magnetic disk through a dedicated DEC MINC 
11/23 Computer. Six groups of 256 samples each 
were taken at each data point. The recorded data 
on the disk were then processed on a DEC PDP 11/70 
minicomputer and plotted on a Benson-Var ian elec- 
trostatic plotter. The available city water did 
not have enough scattering particles to give good 
signal to-noise ratios in the backward scattering 
mode; seeding the test medium with 15.6 pm Dow 
Corning polystyrene particles resulted in signals 
of excellent quality. A TSI hot-film anemometer 
system with a TSI 1231 W conical hot-film probe was 
also used for diagnostic measurements of the jet 
flow at the nozzle exits. A detailed description 
of the experimental apparatus and the flow visua- 
lization and measurement techniques is given in 
Refs. 14-17. 


Test Conditions 


Extensive diagnostic flow visualization studies 
were conducted for equal strength jets and for 
varying heights of the nozzle exit above the ground 
and separation distances between the nozzles. The 
objective was to select a twin-jet impingement 
configuration which has a centrally located, strong 
isolated fountain rising from the ground plate 
without interference from the free jets. Based on 
these flow visualization studies, a normally im- 
pinging, equal-strength twin-jet configuration with 
S/D = 9.0 and H/D = 3.0 was selected for detailed 
LDV measurements. Figure 5 shows the corresponding 
flowfield visualized in a plane connecting the 
nozzle centerlines using the fluorescent dye/laser 
light sheet technique. 


LDV measurements were taken for equal-strength 
jets at twelve heights (Z/D * 0.05 to 2.94) between 
the nozzle exits and the ground. The jet exit 
velocity ^Uj) was 6.71 m/s, resulting in a Re = 

1.70 x 10. All measurements were taken in the 
plane of symmetry connecting the nozzle center- 
lines. The time-dependent streamwise (U + u) and 
cross-stream (V + v) velocities in the X and Z 
directions, respectively, were measured directly by 
the LDV. The mean velocities (U, V) turbulence 


intensities 



Reynolds shear stress 



uv/U 2 ) 


and correlation coefficient 


( — uv/ 



were derived from the basic velocity data. 


Experimental Results 


The experimental data are classified into three 
categories, (a) jet-exit flow; establishing the 
quality of the flow exiting the nozzles (b) jet 
impingement flow (c) fountain upwash flow. 

Jet-Exit Flow Measurements were made in the free- 
jet flow exiting the nozzles to establish unifor- 
mity of the mean velocity profile and the core 
turbulence levels. Both the conical hot-film probe 
and the LDV were used to obtain a cross check of 
the data. Further, the hot-film measurements also 
facilitated spectral analysis of the velocity 
signals . 


1 H9 


Fig. 3 Laser Doppler velocimeter. 




Twin-jet impingement flow with fountain 
formation: S/D ■ 9.0, H/D = 3-0, 

Re = 1.70 x 10. 
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Figure 6 shows the mean-velocity and turbu- 
lence-intensity profiles across the jet operating 
at a Re ■ 1.70 x 10 , taken 1.5 mm below the nozzle 
exit. Here U J is the jet centerline velocity at 
the exit. The jet flow is uniform with a low level 
of turbulence in the core region; the shear layers 
are thin. The agreement between the LDV data and 
hot-film data is good. Because of the large veloc- 
ity gradients across the probe volume, which was 
larger than the shear-layer thickness, the LDV gave 
higher turbulence levels in the shear layers. The 
frequency spectrum of the hot-film signal taken at 
the center of the jet showed no selective or speci- 
fic peaks indicative of unwanted disturbances. 

Jet Impingement Flow The mean-velocity and tur- 
bulence quantities in the jet were normalized with 
Uj. Figure 7 shows the streamwise (U/U ) and 
cross-stream (V/U ) components of the mean veloc- 
ities across the Impinging jet at several stream- 
wise locations, with special emph asis on the stag- 
nation region. The streamwise (V u 2 /UJ and cross- 

stream (1 v /U^) turbulence intensities are shown 

in Figure 8. The Reynolds shear stress (- uv/U 2 ) 
profiles are shown in Figure 9. J 

The measurements, especially the V/U data near 
the ground plate, show a persistent symmetry in the 
flow about the centerline of the jet. The in- 
fluence of the ground plate (jet impingement 
region, Fig. 1) extends to a height of Z/D » 0.75, 
where the jet still has a potential core in the 
streamwise mean velocity profile and the cross- 
stream component of the velocity is close to zero. 
Below Z/D = 0.75, the jet starts deflecting along 
the ground with a rapid decrease in U and increase 
in V. The dip in the U profiles at the center 
seems typical of jet-impingement flows, as also 
shown in the data of Donaldson and Snedeker (Ref. 
18). Figure 8 shows that in the free jet, away 

JT J=f 

from the ground, !u - 21 v ; this trend reverses 


as the jet deflects near the ground. The Reynolds- 
shear-stress profiles (Fig. 9) near the ground show 
substantial regions of zero shear stress (Z/D - 
0.2) despite the large velocity gradients, pos- 
sibly because the turbulence exhibits a delayed 
response to the ground plate, as also evidenced by 

the profiles of and near the ground. 

Limited measurements on the second jet show similar 
initial profiles and jet development. 


fountain upwash 


r iuw 


. fAKui e iu nxustrates tne 

variations of the mean velocities (U and V) across 
the fountain. The distributions of V near the 
ground reveal that the fountain-formation region 
(Fig. 1) extends to Z/D ■ 0.5, where the variation 
In V across the fountain is comparable to that at 
stations farther downstream. However, it should be 
observed that even at Z/D - 0.2 and 0.3, the shape 
of the U profile resembles that typical for other 
downstream stations. The turbulence intensities 

v? and are of the same order of magnitude at a 
given station above the fountain-formation region 
Z/D > °’ 5 ’ In the fountain-formation region Z/D < 

0.5, { 7 was observed to be relatively small 
because of the stagnating flow. The symmetry in 
the data distributions about the centerline of the 


fountain reflects the quality of the present meas- 
urements. The symmetry in the fountain data was a 
major problem in the majority of the earlier inves- 
tigations. 


Similarity of the Fountain Measurements Observa- 
tion of the mean velocity (Fig. 10) and turbulence 
profiles in the fountain suggested self-similarity. 


o,# hot-wire data 
A,0 LDV data 


U/U 



Fig. 6 Mean velocity (U/U ) and turbulence- 

r — ■ — max 

intensity (\u 2 /U ) profiles at the jet 
exit. 
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Variation of (a) streamwise and 

(b ) cr 033-stream mean velocities across 

the impinging jet. 
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The fountain profiles were normalized with U 
and X where U is the local maximum velocity 
and x"?. the haTf-width, is the distance from the 
center' of the fountain to the point where the 
velocity is U /2. U and X were obtained by 
curve-f itting m ?fie fountain mean velocity (U) pro- 
file with a least square curve of the form shown in 
Equation ( 1 ) : 


U 


A + B exp [-(X - 


X ) 2 /2 
o 


s 2 ] 


(1) 


fit 


as suggested by Gilbert (Ref. 4). This curve 
gives the symmetry coordinate X , the maximum 
velocity (A ♦ 3) and the parameter S related to the 
fountain half 

- 1/2 


width X 1/2 by Equation (2). 


1 /2 


2 LOG 


<-)■ 


(2) 


1.177 S (for B»A) 


also 


Fig. 9 
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Variation of the Reynolds shear stress 
across the impinging jet. 


U/U * exp - 0.693 (X/X w .) (for B>>A). 
max r 1/2 

The value of the X for Z/D 2 0.5 was found to be 
very small establishing once again, the inherent 
symmetricity in the present data. Figure 11 shows 
the linear growth of the fountain half-width (X 1/2 ) 
obtained by curve-fitting. The growth rate of 
about 0.16 is lower than that observed in Ref. { ) 
for two-dimensional fountain upwash. Figure 12 
shows the decay of the maximum velocity (U max ) in 
the fountain again obtained by curve fitting. 

Figure 13 shows the streamwise fountain veloc- 
ity (U) profiles shifted to their symmetry point 
and normalized with respect to U and 
obtained by curve fitting. A strflcing similarity 
can be observed in the velocity profiles. Similar 
observations were made for two-dimensional fountain 
in Ref. (4). The profiles below Z/D - 0.5 were 
excluded because they fall into the fountain-forma- 
tion region, although the profiles at Z/D = 0.3 and 
0 2 do not differ much from the similarity form 
shown in Fig. 13. The profiles reach a similar 
form within a short distance (Z/D = 0.5) above the 
fountain-formation region, possibly because of 
greatly enhanced mixing at the base of the 
fountain . 

The variation of the cross stream velocity V 
(Fig. 14) through the fountain was also found to be 
self-similar. These cross-stream mean velocity 
show an expected smooth variation from +ve on one 
side of the fountain to -ve on the other with zero 
crossing at the center of the fountain. The values 
of V at Z/D - 0.5 are slightly higher compared to 
the other downstream stations because of proximity 
to the fountain formation region. 

Figures 15-18 show the various turbulence 
quantities across the fountain nondimensionalized 
with the similarity variables U an * *1/2 ° ~ can 

tained by curve fitting. Once again similarity 
be observed in the turbulence profiles at_various_ 

downstream stations. The profiles of V u 2 and V v 2 
(Figs 15 and 16) show that they are generally of 
the same magnitude and shape. The maximum tur- 
intensities are around 0.5, based on the 
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Fig. 17 Similar profiles of Reynolds shear stress 

(-uv/U 2 ) across the fountain, 

max 



function (-uv/ 
fountain. 




local maximum velocity (U ) , suggesting the 
possibility for even instantaneous flow reversal; 
therefore the LDV is a more reliable Instrument 
than the hot-film anemometer for measuring fountain 
flows. The similarity revealed by the Reynolds 

shear stress (- in?) data (Fig. 17) is significant 
and indicative of the accuracy and precision of the 
experiment because the shear stress data in general 
are particularly sensitive to the measurement 
technique. 

Figure 18 shows the similarity exhibited by the 

correlation function -uv/Vu^ Vv^. The scatter in 
the region beyond X - X 1 / is primarily caused by 
the uncertainty in the measurement of small values 
of the respective turbulence quantities. The 
absolute maximum value of the correlation function 
is observed to be approximately 0.5. 

In general, the turbulence profiles in the 
fountain were observed to take a longer distance to 
attain their self-similar shape than the corre- 
sponding mean velocity profiles. 

Conclusions and Suggested Future Research 

Mean velocity and turbulence profiles were 
obtained at 12 stations across the fountain upwash 
generated by the impingement of two axisymmetric 
jets using the laser Doppler velocimeter (LDV). 

Data were also reported on the impinging jet. The 
turbulence quantities included turbulence inten- 
sities, Reynolds shear stresses and correlation 
functions. The distributions of mean velocity and 
turbulence quantities across the fountain show 
self-similarity when nondimensionali zed with proper 
parameters. The fountain flow was observed to be 
highly turbulent (» 50* turbulence intensity) 
requiring the use of an LDV. The growth of the 
fountain was observed to be linear at a growth rate 
of 0.16. The turbulence and the near ground mean 
velocity measurements in the three-dimensional 
fountain are believed to be the first data reported 
on such flows. 

Additional measurements on the fountain flow 
away from the symmetry plane connecting the nozzle 
centerlines using a three component LDV are sug- 
gested to fully establish the fountain behavior. 

The effect of varying the nozzle separation and the 
height above the ground on the fountain behavior 
need to be established. 
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